Quantitative calculations of charge carrier densities in the depletion layers at 

YBa 2 Cu307_,5 interfaces 
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Charge redistribution at high-T c superconductor interfaces and grain boundaries on the one hand 
is problematic for technological application. On the other hand, it gives rise to a great perspective 
for tailoring the local electronic states. For prototypical (metallic) interfaces, we derive quantitative 
results for the intrinsic doping of the Cu02-planes, i.e. for the deviation of the charge carrier density 
from the bulk value. Our data are based on ab-initio supercell calculations within density functional 
theory. A remarkable hole-underdoping is inherent to the clean interface, almost independent of the 
interface geometry. On the contrary, cation substitution as well as incorporation of electronegative 
impurities can compensate the intrinsic charge transfer and provide access to an exact adjustment 
of the superconductor's doping. The effects of oxygen deficiency are discussed. 

PACS numbers: 73.20.-r,73.40. Jn,74.25. Jb,85.25.Am 
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In recent years, interfaces incorporating a high-T c su- 
perconducting material are subject to a quickly growing 
interest, in both the experimental and theoretical com- 
munity. Due to the possibility of orbital reconstructions 
and a modified chemical bonding at interfaces, design of 
heterostructures with engineered physical properties has 
become a prosperous field. Much effort has focussed on 
YBa 2 Cu 3 07^ based interfaces, often attached to oxide 
compounds [l|, H, Q • In the high-T c cuprates and nick- 
elates, even a very small change of the crystal structure 
or the chemical composition, or the presence of struc- 
tural strain or disorder, can result in fundamental alter- 
ations of the electronic properties [J, [H, [H, 0] ■ From the 
structural point of view, large dielectric constants and 
small carrier densities amplify relaxation effects Q . This 
fact is of special importance for the transport proper- 
ties of the superconductor, since the local charge distri- 
bution in the Cu02-planes is seriously influenced near 
interfaces. Moreover, as a consequence of electrostatic 
screening lengths of a few nanometers and a strong struc- 
tural inhomogeneity, conventional band bending models 
cannot be applied for describing the charge redistribu- 
tion. 

Both the structural relaxation and the local electronic 
states at interfaces of high-T c compounds have been in- 
vestigated for specific situations by various methods, in- 
cluding superlattices composed of different high-T c com- 
pounds [§], doping induced by charge transfer from the 
localized states in a hopping insulator [Io( ]. and cuprates 
in contact to manganites [111 ]. In this context, we have 
recently dealt with the influence of the charge redistri- 
bution at a metal interface on the doping in YBa2Cu3C>7 
(YBCO). It turnes out that the superconducting Cu02- 
planes are subject to a strong intrinsic electron-doping, 
where screening effects arc important for describing the 
charge transfer amplitude [Hj]- The net charge transfer 
amounts to 0.13 electrons per Cu site when the interface 
is oriented parallel to the Cu02-planes, and to 0.09 elec- 
trons for a perpendicular orientation [13| . In addition to 
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FIG. 1: (Color online) Interface between bulk Pd and YBCO 
(precisely: the CuO-chain layer). The Cu02-planes run par- 
allel to the interface. 



this weak dependence on the geometry, the charge trans- 
fer is also robust against the specific metal forming the 
interface, thus giving rise to a "universal" behaviour. In 
this paper we deal with the question, how the intrinsic 
doping changes under different kinds of perturbations of 
the interface. In particular, we give a quantitative evalu- 
ation of various mechnisms that can be used for tailoring 
the charge transfer, which provides a great technological 
potential for material optimization and design. 

While for bulk high-T c compounds a formidable num- 
ber of investigations of impurity effects is found in the 
literature, similar results for surfaces and interfaces are 
rather rare. In addition, studies of interfaces have pre- 
dominately addressed the domain of the superconductor, 
whereas little deal with the interface itself. On the con- 
trary, very much attention has focussed on the influence 
of the O content - in the full range of the phase diagram 
[f3 |. Moreover, it has been observed that the supercon- 
ducting properties of YBCO strongly alter under substi- 
tution of Y by Ca. The critical current at low magnetic 
fields and the critical temperature decrease with the Ca 
concentration (l5j . In the case of grain boundaries, Ca 
doping hence is used to compensate the local charge re- 
distribution [ill. The YBCO crystal structure is subject 
to drastic alterations under the incorporation of defects. 
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FIG. 2: (Color online) CI doped YBCO interface. Every 2nd 
contact metal atom is substituted by CI. The orientation of 
the Cu02-planes is parallel to the interface, which is formed 
by CuO-chains. 
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TABLE I: Selected bond lengths (in A) for different interfaces 
(see the text for details), obtained by structure optimization. 
Bulk YBCO bond lengths are given for comparison. 



Especially, the lattice parameters change considerably as 
a function of the O stochiometry [T3|. Experiments in 
the late eighties had already indicated that the insertion 
of F/Cl into underdoped YBCO can enhance the critical 
temperature where the F/Cl atoms tend to occupy 
the empty metal sites in the CuO-chains [HI, [2(| ■ Inter- 
estingly, this insertion of halogenides is accompanied by 
a higher hole content in the Cu02-planes, which suggests 
that a F/Cl impurity has the potential to re-extract the 
excess charge arising at an YBCO interface. 

Our data are obtained from density functional theory, 
applying the generalized gradient approximation and the 
Pcrdew-Burkc-Ernzernhof scheme. We use the WIEN2fc 
program package [2ll ] , with a mixed linearized augmented 
plane wave and augmented plane wave plus local orbital 
basis. Moreover, we set RK maK — 7 when we investigate 
CI impurities and RK max = 4.7 otherwise. The radii of 
the muffin tin spheres are (in Bohr radii): 1.8 for Cu, 2.2 
for CI, 2.25 for Ca, 2.3 for Y, 2.5 for Ba, 1.45 to 1.55 for 
O, and 2.1 to 2.4 for Pd. 

Because band-bending takes place on the length scale 
of the YBCO c lattice constant, the electronic structure 
of an interface is accessible to a supercell approach with 
periodic boundary conditions [13] • It is convenient to use 
as the metallic substituent fee Pd, which has a minimal 
lattice mismatch of only 0.7% with respect to the YBCO 
a/6 lattice constant. We have checked that our results do 
not depend on this choice by comparing the charge trans- 
fer at YBCO-Pd and YBCO-Ag interfaces, for which we 
find almost identical data in all relevant respects [H, E3] • 
For illustrating the principle of the supercell setup, Fig. 
[1] displays the near- interface region of the supercell for a 
clean YBCO-Pd interface oriented parallel to the Cu02- 
planes. It comprises two YBCO unit cells terminate d by 
the CuO-chain layer, as indicated by experiment [23L[24|. 
and four Pd fee unit cells stacked in [001]-direction. The 
supercells under consideration all have a similar size. 

We use the experimental YBCO bulk lattice constants 
a = 3.865 A and b = 3.879 A (25|, which is likewise ap- 
plied to the Pd region in order to avoid lattice mismatch 
in the contact plane. In the perpendicular direction, the 
standard Pd lattice constant c = 3.89 A is an adequate 



choice. By the minor dependence of the charge transfer 
on the interface geometry, all subsequent results refer to 
the setup of Fig. [TJ In the case of fully oxydized YBCO, 
one half of the topmost Pd atoms consequently has an O 
bonding partner in the interface CuO-chain layer, while 
the other half is located next to an unoccupied site. We 
count atomic layers with respect to the interface, i.e. the 
first Pd layer is attached to the first atomic layer of the 
YBCO cell, formed by the the CuO-chains. 

The following 3 mechanisms of tailoring the intrinsic 
doping are investigated: First, the incorporation of elec- 
tronegative ions is realized by substituting each second 
interface Pd atom without an O neighbour by a CI atom. 
This arrangement is shown in Fig. [2j after structure op- 
timization. Second, cation substitution is studied for the 
replacement of Y by Ca, where we choose a Ca content 
close to the optimal value of w 30% for grain boundary 
doping [16j. To keep the computational effort manage- 
able, we replace each fourth Y atom by Ca, i.e. we deal 
with Yo.75Cao.25BCO. Third, the effects of O deficiency 
are addressed for YBCO6.75, where each fourth interface 
O atom, all from CuO-chains, is removed. We note that 
the system therefore stays far in the metallic regime [26j j , 
so that valid results are expected from a band structure 
calculation using the present approximations. 

Turning to the chemical bonding in the vicinity of the 
interface, there is always the same, expected behaviour. 
We find a strong tendency towards Pd-0 bonding and a 
Cu-Pd repulsion. The Cu-Cl repulsion is weaker. More- 
over, the CI incorporation induces only weak structural 
distortions of the superconductor, whereas the first two 
Pd layers are subject to a finite relaxation. In fact, the 
relaxation in the superconductor domain declines quite 
quickly for all interfaces under investigation, so that the 
second Cu02-plane is virtually undistorted. This is par- 
ticularly true for the O deficient case where we may an- 
ticipate a different behaviour. Selected bond lengths are 
compared to each other in Table [TJ Reflecting covalent 
Pd-0 bonding across the interface, the nearby Pd sites 
are subject to shifts of the valence electronic states from 
lower energy to about — 2eV, with respect to the Fermi 
energy. This gives rise to a pronounced structure in the 
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FIG. 3: (Color online) Partial Cu 3d DOS for the 2nd Cu0 2 - 
plane off the clean (upper panel) and CI doped (lower panel) 
YBCO interface, as compared to the bulk DOS. 



Pd Ad density of states (DOS), characterized by a signif- 
icant Pd-0 hybridization. However, even for the second 
Pd layer, the DOS resembles the bulk Pd 4d DOS very 
well, and a nearly perfect agreement is obtained for the 
third layer. As a consequence, structural distortions are 
screened quickly on both sides of the interface, which a 
posteriori justifies the size of our supercells. 

We first investigate the effect of the CI doping on the 
YBCO electronic states. Cl-Pd hybridization is tiny for 
atoms in the first Pd layer, according to an ionic bond- 
ing [27], [H[. However, there are remarkable admixtures 
of the nearest neighbour Pd atom (see Fig. (2) in the CI 
dominated DOS range, which we attribute to enhanced 
bonding because of a very short distance. Of course, it is 
to be expected that the CI extracts charge from the sur- 
rounding area until a CI - ion state is reached. Since no 
reduction of the Pd charge is observed under CI doping, 
the extraction should affect only the YBCO region. As- 
suming that the charge reduction affects all 7 Cu atoms 
per unit cell to a similar degree, naive electron counting 
indicates that the intrinsic charge transfer due to the in- 
terface (0.13 electrons per Cu site) can be compensated 
by this mechanism quite well. 

We therefore compare in Fig. [3] the DOS for the clean 
and for the CI doped interface to bulk YBCO. The data 
refer to the second Cu02-plane, because here the struc- 



FIG. 4: (Color online) Partial Cu 3d DOS for the 2nd Cu0 2 - 
plane off the Yo.7sCao.25BCO (upper panel) and O-deficient 
YBCO6.75 (lower panel) interface. In each case, the findings 
are compared to the corresponding DOS of bulk YBCO. 



tural relaxation has already decayed. According to Fig. 
[21 the Cu 3c? DOS confirms our speculation. Whereas a 
strong reduction of the (unoccupied) states is visible for 
the clean interface, this charge transfer is fully compen- 
sated for the CI doped case. The areas under the DOS 
curves in the lower panel of Fig. [3] coincide almost per- 
fectly, reflecting identically doped Cu02-planes. In fact, 
the suppression of charge transfer to the YBCO domain 
due to the CI atoms is surprisingly strong. Electronega- 
tive impurities consequently seem to form a rather rigid 
potential barrier which prohibits any intrinsic doping. 

Another way for overcoming the hole-underdoping at 
interfaces is simultaneous overdoping: replacing Y 3+ by 
Ca 2+ . The initial doping state should be restored, since 
the introduction of Ca increases the hole concentration. 
The upper panel of Fig. 2] confronts the DOS of the pure 
YBCO interface with the results for an Yo.7sCao.25BCO 
interface. Whereas differences in the Cu 3d DOS shape 
reflect remarkable alterations of the electronic states, in 
particular of the Cu-0 bonding, their occupation fulfills 
exactly the trends mentioned before. The interface dop- 
ing is even overcompensated for the configuration under 
investigation, which is visible in Fig. 0] by a shift of the 
DOS curve to higher energies. Further evaluation of the 
data gives a charge reduction of 0.08 electrons. 25% Ca 
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doping thus comes along with an injection of 0.21 addi- 
tional holes per Cu site into the Cu0 2 -planes. 

The fact that local Ca doping of grain boundaries can 
enhance the supercurrent density has been attributed to 
a reduction of the charge accumulation at grain bound- 
ary dislocations, which causes a hole depletion [29|. The 
size of the perturbed regions seems to be smaller for Ca 
doped samples and the interface conductivity thus is im- 
proved. Resistivity measurements in fact reveal changes 
of the height and the sha pe o f the electrostatic potential 
barrier under Ca doping [301 ] - However, the experiments 
likewise indicate that such modifications cannot explain 
the observed effect of Ca doping on the transport prop- 
erties. It therefore has been speculated that extra holes 
must be generated for some reason, which finally is sub- 
stantiated by our findings. A grain boundary appears to 
be quite well described in terms of an interface between 
an unperturbed superconductor domain and the charge 
reservoir of the grain boundary dislocations. Ca doping 
then compensates the intrinsic charge transfer present at 
this interface and so introduces additional holes. 

As is known from many experiments, O deficiency has 
a very strong influence on the hole content of the Cu02- 
planes, see [311 ] . for instance, and the references therein. 
In analogy to our previous proceeding, we now quantify 
this effect. Evaluating the DOS given in the lower panel 
of Fig. 21 we obtain a reduction of the hole count in our 
YBCO6.75 interface by 0.10 holes per Cu site. However, 
this reduction is not due to charge transfer, as we do not 



observe an energetical shift of the DOS curve, but traces 
back to a modified Cu-0 bonding. 

We have discussed electronic structure calculations for 
prototypical interfaces of a high-T c superconductor and 
a normal metal. In each case, we have evaluated quanti- 
tatively the interface charge transfer from the metal into 
the YBa2Cu307-5 domain. The intrinsic underdoping of 
0.13 holes per Cu atom (for 5 = 0), is fully suppressed 
by electronegative impurities. Cation substitution in the 
superconductor likewise leads to a (over-)compensation, 
where the initial hole concentration is restored for a 15% 
Ca content. Since O deficiency at the interface has just 
the opposite effect, it is therefore possible to exactly de- 
sign the doping state of the superconductor. Reduction 
of the O content by S = 0.25 leads to an underdoping of 
0.10 holes per Cu atom. We expect that all these values 
likewise apply to the intentional doping of grain bound- 
aries, which paves the way for a systematic optimization 
of transport properties. 
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